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Abstract 

A model-independent analysis for the exclusive, rare B — > ir + ir~l + l~ decay is 
presented. The dependence of the various physically measurable asymmetries and 
CP violating asymmetry on the new Wilson coefficients is studied in detail. It is 
observed that different measurable quantiles are very sensitive to the different new 
Wilson coefficients, i.e., to the existence of new physics beyond standard model. 
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1 Introduction 



Rare B meson decays, induced by flavor-changing neutral current (FCNC) b — > s or b — > d 
transitions, occur at loop level in the standard model (SM) and is a potential precision 
testing ground for the SM, as well as attracting the theoretical interest as a promising tool 
for establishing new physics beyond SM. New physics effects show themselves in rare B 
meson decays in two different ways, namely, through the Wilson coefficients which could 
be distinctly different from their SM counterparts or through the new structures in the 
effective Hamiltonian (see for example [01 )• One of the main goal of the working 
B factories and of the future hadron colliders is the study of the rare B meson decays. 
Therefore theoretical and experimental investigation of the rare decays of B mesons receive 
special attention. The aim of this paper to investigate the rare decays B — > K~n + £ + £~ and 
B — > tt + tt~£ + £~ when Kir and tttt systems are the decay products of K* and p mesons, in 
a model-independent manner. 

The inclusive b — > s£ + £~ and exclusive B — > K*£ + £~ decays were analyzed in a model- 



independent way in [I ![ and [12| , respectively. The present work is an extension of the 
previous studies of the B — > K*(p)£ + £~ decay. The distribution in cos^p, where 6 P is the 
angle of the K~(tc~) in the K~7t + (tt~7t + ) center of mass frame, and the dependence on 
the azimuthal angle if between the £ + £~ and K~7t + or ir~7t + planes, which does not exist 
in the B — > K*(p)£ + £~ decay, can provide additional new information. This information 
is sensitive to the polarization state of the vector meson K*(p) and therefore opens new 
possibility in probing the structure of the effective Hamiltonian. Angular distributions and 
CP asymmetries in the B — > K~n + £ + £~ and B — > tt + it~£ + £~ decays in the SM were studied 



m 



The paper is organized as follows. In section 2, we present the general, model indepen- 
dent expression of the decay distribution for the B — > Kh^-k)^ £~ decay, in terms of the 
helicity amplitudes, including the non-zero lepton mass effects. In section 3, we present 
our numerical analysis for different angular distribution together with a brief concluding 
remark of our results. 



2 Theoretical background 

The matrix element for the B K*(p)£ + £~ decay at quark level is described by the 
b — > f£ + £~~ (f = s,d) transition. Following the works [|T^]-[|T^], the matrix elements of the 
b — > f£ + £~ decay can be written as the sum of the SM and new physics contributions as 



= V tb V t * f | Csl /«v |r L b l lv !L + C BR fia, v q - R b l^i 

+ cti hi,b L l Ll n L + ct l R hi^L titf in + c RL f Rl ,b R l Ll n L 

+ C RR fn'y^bn Ir^Ir + C LRLR f L b r £l£r + Crllr Jr^l £l£r 

+ C LRRL f L b R £ R £ L + Crlrl f R b L IrIl + C T f^b la^l 

+ iC TE e^fa^b£a a A , (1) 
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the first four terms contain the SM contribution too, while all the others describe the new 
physics effects and where 



75 /I + 75 



L(R) 



In further analysis we will neglect the tensor type interactions (terms ~ CV and Cte since 
physically measurable quantities are not sensitive to the presence of tensor type interactions, 
as is shown in |I| [Tl| ) . 

The matrix elements for the precess B — > PP'£ + £~ can be obtained from the matrix 



element B — > V£ + £ (V = K*, p) in the following way 

M = M B ^ V e+i- Tl(s M )(n + (p + )n-(p_)\V( Pv , A)) 



(2) 



where we assume that the resonance contribution of the intermediate vector meson can be 
implemented by the Breit-Wigner form 



n(s M ) 



%- m v + im v T\ 



(3) 



where % = {p + +pJ) 2 and m v and T v are the mass and width of the vector meson (K*, p). 



For the decay part of the vector meson we use |15 



(P(p + )P'(p.)\V(p v ,X)} = VBRY x x : je P ^) 



(4) 



where Y[ n {9p ) (p) are the J = £ spherical harmonics, and the angles 9p and tp belong to 
those of the final P meson in the vector meson's rest frame. The coupling of V — > PP' 
decay is effectively taken into account by the branching ratios. 

From (HP we see that in order to calculate the matrix element B — > PP'£ + £~ , the matrix 
element of the B — > V£ + £~ decay is needed in the first hand, for which the matrix elements 



V 



/7 M (1± 75 )6 B{p B \ 



V 



/«V(l + 75 )& B(p B )) and (V f{l±j s )b B(p B )) , (5) 



need to be calculated, whose most general forms can be written as 
[V(p v ,e) / 7Ai (l±7 5 )6 B{p B )) 



j. 2V (q 2 



m B + my 



± ie* (m B + m v )A 1 (q 2 ) =f i(p B + p v )^*q) 



m B + my 



2m 



T«V 2 



v 



e*q) A 3 (q 2 ) - A (q 2 



(6) 



V(p v ,e) fia fiu q v (l± l5 )b B{p B ) 



^ X(7 e*"p v q' 7 T 1 (q 2 )±2 l 



±2i(e*q) 



q 2 



s M ) - (PB + Pv)^( £ *q) T 2{q 2 ) 



q»- {pb +pv)ti 2 

m% - s M 



Uq 2 



(7) 



where q = p B — pv is the momentum transfer and e is the polarization vector of the vector 
meson. In order to ensure finiteness of (§) at q 2 = 0, we demand that A 3 (q 2 = 0) = A (q 2 = 



2 



0). The matrix element (V /(l ± 75)6 Bj can easily be calculated from Eq. (|6]). For this 
aim it is enough to contract both sides of Eq. with q p and use equation of motion. 
Taking mass of the strange quark to be zero, it gives 



V(pv,e) /(1± 75 )6 B(p B )\ = —{+i(e*q)(m B +m v )A 1 



± i(m 2 B - s M )(e*q)- 



A,{a 2 ) 



±2tm v (e*q)[A 3 (q 2 ) - A (q 2 )} } . (8) 



m B + my 

Furthermore, using the equation of motion, the form factor A3 can be expressed in terms 
of the form factors A\ and A2 (see ||16|| ) 



A 3 (q 2 



m B + m v 



A x { 



m B — my 



A,(a 2 ) 



(9) 



2m v 2my 
Using this relation, the matrix element (§) can be written in the following form 

(V(py,e) f(l± l5 )b B(p B )) = —{T2imv(e*q)A (q 2 )} . (10) 
x ' m\j L J 

As a result of the above considerations, the matrix element of the B —>■ V£ + £~ decay can 
be determined straightforwardly 

Ga 



M(B -> Vt 



—r-^ti 



x 1^(1 - <y 5 )£ [ - 2V Ll ^ vpa e* v p v q a - tV L2 e; + iV L3 (e*q)(p B + p v ) p + iV LA {e*q)q» 

+1^(1 + l5 )£ [ - 2V Rl e flupa e* u p p v q CT - iV^e* + iV R3 {e*q){p B + p v ) p + iV RA (e*q)q p 

+1(1 - l5 )£[iS L (e*q)] + 1(1 + l5 )£[iS R (e*qj\ j , (11) 

where Vx 4 and V Ri are the coefficients of left- and right-handed leptonic currents with vector 
structure, and Sl,r are the weights of scalar leptonic currents with respective chirality, 
respectively, whose explicit forms are given as 



Vl 4 



(C L °l + C RL ) 2(C BR + C S L)-k , 

m B + m v q 2 

(Cl°l - C RL )(m B + m v )A 1 - 2(C BR - C SL )^(m 2 B - s M 



ClL ~ CRL -A 2 - 2(C BR - C SL ) 1 



m B + my 



r 



T 2 + 



q 



(C L °l - C RL )^f(A 3 - Ao) + 2(C BR - C SL 



m B - s M 

2 ' 





= v Ll (c?i- 


~* b LR > 


C R L — 


c RR ) 




= vuctl- 


/~itot 
b LR 5 


C R L — 


C RR ) 




= v L ,(ctt- 


s~itot 
°LK 5 


Crl — 


c RR ) 




= V U {C%- 


sitot 


C R l — 


c RR ) 


S l 


= —(Clrrl 


— C R l R L 




A)) , 






V m 6 
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A' 



— (Clrlr ~ C R llr) ^ V A[) 



m 6 
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In order to obtain the full helicity amplitude of the B — > PP'£ + £ which follows from Eq. 
(§), the helicity amplitude of the B — > V£ + £~ decay must be written, which we denote as 



KA^t^t — \" , n , jX 1 X 1 ttX 1 



(12) 



where 



■ X(\i 

V* 

jjXiXi 

Xy* 



(vipv^^^BipB: 



J, 



(13) 



where £y* is the polarization vector of the virtual intermediate vector boson satisfying the 
relation 
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EfJL *V 



where the summation is over the helicities Ay* = ±1, 0, s of the virtual intermediate vector 
meson, with the metric defined as r}± = rjo = —r] s = 1 (for more detail see |L7|, |18|). In Eq. 
(0), and J 1 are the leptonic and hadronic currents, respectively. Using Eqs. (]TT|) — (]T3|) 

for the helicity amplitudes M.\. 1 we get the following expressions 



M ± - 

Mt 
Ml 

-Mn" 



sin#A+ + , 
(-l±cos0)A^ 
(1 ±cos6)A ± + . 
sin6L47~ , 



B 



cos 6A£ 
sw.6A£~ , 
sin6'yl + , 
cos6Aq~ + B~ 



(14) 
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where 



H ± + 



C RL (l-v) + C RR (l + v) 



h 



± 



A+ + = ±V2m t {(C% + Cf R )E ± + - (C BR G ± + CW) + (C RR + C RL )/i ± } , 
A ± = — A± + , 

Al- = ^{[Cl°i(l-v) + Ci° R (l + v) 

+ — (C BR G± + C SL g±) } , 
A~ ± + = A + ± -(v^-v) , 

A++ = 2m e [(C™ + C? R )H Q + (C RL + C RR )h + ^ (Crf + C S l9o) 

Aq — — Aq + , 

5 ++ = -1m s {[Cf R - Cii)H° s + (C RR - C RL )h° s } 

- —q 2 [(1 - v){C LRLR - C RLLR ) - (1 + v){C LRRL - C R l R l ) J ii/g | , 
fl ~ = - -v) , 



(15) 



" = [ClK 1 ~ v ) + Crr(\ + V )] H ° + [C M (1 - u) + Chh(1 + v)]h 

2 

+ ^ (CbrGq + Csl9o) } , 

where superscripts denote helicities of the lepton and antilepton and subscripts correspond 
to the helicity of the vector meson (in our case p or K* meson), respectively. 

H ± = ±\ 1 /\m 2 B , s M , q 2 ) V{q2) + (m B + m v )A 1 {q 2 ) , 

rriB + my 

I [ 

- (m| - s M - q 2 )(m B + m v )A 1 (q 2 ) 



H, 



2y/s M q 2 

w 2 2 ^ Mq 2 ) 

\{m B ,s M ,q 



tjo X 1/2 (m 2 B ,s M ,q 



2 



+ 2^(A 3 - A ) 
G± = -2 
G 



j\ ( ^ 1 

- (m B + m v )A 1 (q 2 ) H (m B - s M ) 

m B + my 



(16) 



1 



\fs M q 2 
q 2 



± A 1/2 (m|, s M , g 2 )T!(g 2 ) + (m| - s M )T 2 (g 2 ) , 

m B - %)K - s M - q 2 )T 2 {q 2 ) - X(m 2 B ,s M ,q 2 )(T 2 (q 2 ) 

Uq 2 ) 



+ 

m B - s M 
h± = H±(A 1 -> -Ai, A 2 -> -A 2 ) , 
h° s = HKA! - -Ax, A 2 - -A 2 ) , 
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where is the polar angle of positron in the rest frame of the intermediate boson with 
respect to its helicity axis. Note that we take py = sm, but not m v , in order to take into 
account Vs being a virtual particle which subsequently decays into 7r + 7r~ or K~n + pair. 
Remembering that the existing CLEO result [|IIJ for the P — > X s ~f and P — > K*^y decays 
impose strong constraint on the parameter space Cbr and Csl- For this reason here we 
assume that they are equal to each other in the SM. Hence we will take 



LL 
U LR 



r< e ff 



— ClO + C-LL i 

9 + Cio + Clr ■ 

Using the expressions of the helicity amplitudes for the differential decay rate width of the 
B -> K(-> PP')£ + £- decay we get 

3G 2 a 2 2 

= — 3 o VtbV t *e cIsm dq 2 d cos 9 p d cos 9 dtp 

2 i '7T b m B s M q 

x A 1/2 (m|, s M , g 2 )A 1/2 (s M , m|, m|,)A 1/2 (g 2 , m 2 , m 2 ) 



myry/7T 



■B(V -> PP') 



x ^ 2 cos 2 P cos 2 JVi + sin 2 iV 2 + 2 cos Re[iV 3 ] + iV 4 



+ sin 2 P sin 2 9 N 5 + ( 1 + cos 2 0) iV 6 + 2 cos iV 7 + 2 sin(2<^) sin 2 9 lm[N 8 



- 2 cos(2<^) sin 2 Re[iV 8 ] + V2 sin(20 P ) sin 6 1 cos ip Re[ cos 9 Ng + N- 



(17) 



where 



iV 2 
iV 3 
iV 4 

iV 5 
iV 6 

iV 7 

iV 8 
iV 9 



i(20 P ) sin 


sin p Im[ cos N n + N 12 ] j 


AT 


2 




2 


+ 




j 




2 




2 




+ 


A + 


5 



Pn 



AT" P 



AT 



+ 



+ 



AT + 
[At 
A+ + (At 
+ A 

+ B 
N n = N 9 (A 
N 12 = N 10 (A 



B 

At 
AZ 
A~ 



2 

+ 




2 
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2 
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A + _- 


2 

+ 
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(At 


I*- 


f A- 



P +H 

-°o 



AZ~ 

++ 
+ 



■A? 

K~ 

■A? 
-A 



- A + - 



+ At 



A' 



A_ 
At 



+ a: 



- A-+ 



A_ 
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+ AS 
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A: 
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.4 
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-A] 



;is) 
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where Op is the polar angle of the pseudoscalar P meson momentum in the rest frame of the 
vector meson, with respect to the helicity axis, i.e., the outgoing direction of V meson, and 
(p is the azimuthal angle between the planes of the two decays V — ► PP' and V* — ► £ + £~ . 
Kinematically allowed region of the variables are given as 



[mp + m P i) 2 < Sm < ( m B — 2m^) 2 

Am 2 — ? 2 — ( m B ~ 
-1 < cosO P < 1 
-1 < cos6 < 1 , 
< ^ < 2tt . 



(19) 



We note that, in further analysis the narrow-width approximation for V meson will be 
used, i.e., 



lim 



T v m v /iT 



v^o (s M ~ m 2 v ) 2 + m\Y\ 



S(s M ~ m 



2 ) 
v) > 



by means of which integration of Eq. fll7|) over sm can easily be carried and the differential 
decay rate with respect to dilepton mass q 2 , azimuthal angle (p, polar angles 6p and 6 can 
be written as 

3G 2 a 2 



dT 



B{V -> PP') dq 2 d cos Op d cos d<p 



2 17 7r 6 m Bj m 2 / q 2 

\l/2/ 2 2 2\\l/2/ 2 2 2 \\l/2/ 2 2 2\ 

x A 1 {m B , my, q )X ' (m v , m P , m P ,)\ 1 (q , m|, rrq) 

x 1 2 cos 2 0p [ cos 2 9 Ni + sin 2 6 N 2 + 2 cos 6 Re [iV 3 ] + iV 4 

+ sin 2 6 P [ sin 2 6 N 5 + (l + cos 2 0) iV 6 + 2 cos iV 7 + 2 sin(2<^) sin 2 Im[A^ 8 
- 2cos(2(^)sin 2 0Re[iV 8 ] + V2sm(26 P ) sin cos ip Re [ cos 9 Ng + Niq] 



(20) 



- V2sm(26 P ) sin sin cp Im[ cos # iVn + iV : 
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for which we will use the experimental results for the V — > PP' namely, i3(p — > 7r + 7r _ ) = 
B(K* — > i^7r) = 1. It should be noted here that in addition to the variables that exist in 
B — > V r £ + £~ decay, there appears a new variable Op, and contrary to the B — > case, 
the dependence of the cascade decay B — > V"(^ PP')£ + £~ on the azimuthal angle <£> is 
not trivial. Incidentally, we would like to remind the reader that, if Eq. (|20D is integrated 
over 6p and <p, the differential decay rate for the B — > V r £ + £ _ decay is obtained. The 
model independent analysis of the B — > K*£ + £~ decay is presented in [12], in which the 



dependence of the experimentally measured quantities, such as branching ratio, forward- 
backward asymmetry and longitudinal polarization of the final lepton and the ratio T^/IY 
of the decay widths when K* meson is longitudinally and transversally polarized, on the 
new Wilson coefficients, are systematically studied. As has been noted already, our main 
goal in this work is to investigate the dependence of such angular distributions on the 
new Wilson coefficients in the B — > V(V — > PP')£ + £~ decay, which does not exist in the 
B — > V£ + £~ decay and not studied in ||12|| . 

It can easily be seen from Eq. (p0|) that the cascade decay B — ► V(— > PP')£ + £~ has a 
rich angular structure. Therefore, in the light of this observation, a thorough investigation 
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of different distributions will prove useful in separating various angular coefficients experi- 
mentally. Along the lines as suggested by [fOfl , we adopt two different strategies in further 
analysis of the problem under consideration, namely investigation of the various individual 
angular distributions and asymmetries and their relation to the new Wilson coefficients. 

For the purpose of studying single angle distributions, we integrate Eq. fl20|) over q 2 , 9 
and (p, which takes the form 



dT 



d cos 9p 



2 cos 2 6 P 



l Nl + ^N 2 + 2 N 4 



+ sin 2 9 P 



- 3 N i + + -N e 



(21) 



where we introduce the notation Ni = J Nidq 2 . Defining an asymmetry parameter ag p , 
from the angular distribution W^cos^p) = f + ag p cos 2 dp we get 



OiQ, 



N\ + 2N 2 + 3A^ 4 



f . 



N 5 + 2N 6 

Integrating Eq. (^) over q 2 , dp and ip, for the polar cos 9 distribution we get 
dT 4 



(22) 



d cos 9 3 
4 

+ 3 



cos 2 9 ^ + sin 2 9 N 2 + 2 cos 9 Re[N 3 ] + N 4 
sin 2 9 % + (f + cos 2 6) N 6 + 2 cos 9 N 7 



from which one can write the angular distribution W{cos9) = 1 + otg cos9 + (3$ cos 2 9, with 
the asymmetry parameters ag and (3g being defined as 



06 



2Re[iV 3 ] + 2N 7 
N 2 + N 4 + N 5 + N e 
% - N 2 - N 5 + N 6 
N 2 + N 4 + N 5 + N & 



(23) 



Finally we consider the azimuthal angle ip distribution, which is obtained by integrating 
Eq. fl20D over the parameters q 2 , 9 P and 9 to yield 



dT 4 

dp 3 
4 

+ 3 



l Nl + ^N 2 + 2N 4 



4 



- N 5 + - N 6 + - sm(2<p) hn[N g ] - - cos(2<^) Re[iV 9 ] 



and the azimuthal angle p asymmetry parameters and a v can be extracted from W(<p) 
1 + a v sm{2p>) + /3 V cos(2(^) to give 



0u 



4lm[N 9 



iVi + 2N 2 + 3N 4 + 2N 5 + 4N 6 

^ -4Re[iV 9 ] ^ 
Ni + 2N 2 + 3N A + 2N 5 + AN 6 



(24) 
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A second strategy for separating various angular coefficients experimentally, is to define 
suitable asymmetry ratios that project out the partial rates from Eq. (|20|). For this purpose 
we consider the following asymmetries (see also |T8| ) 



where 



dl» - dT(<p + tt/2) + dT(<p + tt) - dT(<p + 3tt/2) 
v ' dT{<p) + dT( V + tt/2) + cir(^ + tt) + dr(yj + 3vr/2) ' 1 j 

7T TT 
< tn < _ 

4 - ^ - 4 ' 

M = (26) 



AT = dr{6,6p,(p)-dr{6,6p,ip + 'K)-dT{6,iT-6p,ip)+dr(6,iT-8p,ip + 'K) 

- dT(ji - 9, e P , <p) + dT(ir - 9, e P , <p + tt) + rfr(vr - 9, tt - P , y>) 

— rfr(7T — 9,TT — 9p, (p + 7r) , 



< P < TT/2 , 

7T 7T 
< w < 

2 - ^ - 2 ' 

and the denominator D is given by the same expression, with plus signs everywhere, 

rfr(g P , g) - rfr(0 P , <p + tt) - rfr(7r - gp, y) + dr(7r - 9 P , <p + tt) 

< 9 P < tt/2 , 
-| < V? < vr/2 . 



In expressions (|25|)-(p7|) the angles that do not appear in the arguments of the differential 
rate dT have been integrated out over their physical ranges (0 < 9, 9 P < tt, < ip < 2tt). 
Integrating over all variables, we are left with the expressions A^, Ai, A2, which depend 
only on the Wilson coefficients, as follows (here ~ in the notation refers to integration's 
being performed over all variables) 

A = _ o-iic[i v 9 j _ ^ , 2i 

V tt(n 1 + 2N 2 + 3N 4 + 2N 5 + 4N 6 ) ' 

Ai = ^m^L ^ , (29) 



A 2 = „ „ — . (30) 

y/2 [Nt + 2N 2 + 3iV 4 + 2N 5 + 4AV 



-8Re[A^ 9 ] 


tt (jVi + 2N 2 + 3N 4 + 2N 5 


+ 4iV 6 ) ' 


-V2Re[N w ] 




tt (jVi + 2N 2 + 3iV 4 + 2N 5 


+ 4iV 6 ) ' 


3Re{N n ) 





Before proceeding further, we would like to consider the CP-violating observables that can 
be constructed by combining the information on B and B decays, namely, B — > PP'£ + £~ 
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and B — > PP'£ + £ , and define CP-odd asymmetry in the following way 

A CP = ^£ , (31) 

where V and T are the decay widths of the B — > PP'£ + £~ and B — > PP'£ + £~ processes, 
respectively. Explicit form of V can easily be obtained from Eq. (^), by performing 
integration over the variables q 2 , 6, 8p, and <p. The decay width for the conjugate process 
can again be obtained from Eq. ( p0|) by making the replacement iVj — > iVj, where iVj are 
the functions for the conjugate processes. It should be noted that we consider a case in 
which all new Wilson coefficients are real. Furthermore, form factors which enter into Eq. 
(|I~6|) are computed in framework of light cone QCD sum rules method |[22| and this 



nonperturbative approach predicts that all form factors are real as well. In other words, 
there is no any new source for CP violation other than that are present in the SM. As is 
well known, for the B — > p£ + £~ decay in the SM only the coefficient Cg contains both 
a weak phase tp w (associated with the imaginary part of the CKM matrix element) and a 
strong phase 5s (attributed to the imaginary parts of the cc and bb loops). As a result of 
these considerations, it follows then that the decay width for the conjugate B — > PP'£ + £~ 
process can be obtained from B — > PP'£ + £~ channel by the replacements ipw —> —<£w an d 
5 S -> Ss- 



3 Numerical analysis 

In this section we present our numerical results for the asymmetries Aqp, ag p , ag, j3g, 
(3^ , A v , a<p , Ax , A<p and A 2 for the exclusive rare B — > -n + n~£ + £~ decay only. We take 
hadronic form factors from Table I and the Wilson coefficients from Table II. The values of 
the main input parameters used in our analysis are: = 4.8 GeV, m c = 1.35 GeV, m p = 
0.77 GeV, m T = 1.78 GeV, = 0.105 GeV and m B = 5.28 GeV. Here we note that 
the results for B — > Ktt£ + £~ can be obtained from B — > tttt£ + £~ by replacing Vtd with Vt s , 
B — > p transition form factors with B — > K transition form factors, and with mx- 

In the present work we choose light cone QCD sum rules method predictions for the 
form factors. In our numerical analysis we will use the results of the work [J2T], in which 
the form factors are described by a three-parameter fit where the radiative corrections 
up to leading twist contribution and SU(3)-breaking effects are taken into account. The 
g 2 -dependence of the form factors, which appears in our analysis could be parametrized as 



1 — a F s + bp s 2 



where s = q 2 /m 2 B is the dilepton invariant mass in units of B meson mass, and the param- 
eters F{0), a F and b F are listed in Table 1 for each form factor. In the SM the Wilson 
coefficients C7 ~ {nib) and Cio(mb), whose analytical expressions are given in [23, 24 



arc 



strictly real as can be read off from Table 2. In the leading logarithmic approximation, at 
the scale 0(/i = m b ), we have 

C e 7 ff {m b ) = -0.313, 
Ctl f {m b ) = -4.669 . 
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F(0) 


dp 


b F 






0.372 


1.40 


0.437 


aT 




0.261 


0.29 


-0.415 


At 


->p 


0.223 


0.93 


-0.092 


yB- 


^p 


0.338 


1.37 


0.315 


T B- 


^p 


0.143 


1.41 


0.361 


Tt 


->p 


0.143 


0.28 


-0.500 




->p 


0.101 


1.06 


-0.076 



Table 1: The form factors for B — > pvl in a three-parameter fit 



Ci 


c 2 


c 3 


c A 


c 5 


c 6 


c 7 e// 


c 9 


r eff 

°10 


C ( ) 


-0.248 


1.107 


0.011 


-0.026 


0.007 


-0.031 


-0.313 


4.344 


-4.669 


0.362 



Table 2: The numerical values of the Wilson coefficients at \i ~ m b scale within the SM. 

Although individual Wilson coefficients at /i ~ mj, level are all real (see Table 2), the 
effective Wilson coefficient Cg (mj, s) has a finite phase, and in next-to-leading order 



Cg ff (m b ,s) = C 9 (m b ) 



a 8 (fi) 
1 H wis 



7T 



+ Y SD (m b , s) + Y LD (m B , s) 



(32) 



where C$(m b ) = 4.344. Here u (s) represents the 0(a s ) corrections coming from one-gluon 
exchange in the matrix element of the corresponding operator, whose explicit form can be 
found in J23J. In (|3*2]) Y SD and Y LD represent, respectively, the short- and long-distance 
contributions of the four-quark operators O i= i t ... t Q [^3], p4fl . Here can be obtained by a 
perturbative calculation 



Y SD (m b ,s) = g(m c ,s)C®--g(l,s)[4C 3 + 4C i + 3C 5 + C 6 \ 

- ^g(0,s)[C 3 + 3C 4 ] + ^[3C 3 + C 4 + 3C 5 + C e ] 

- A„[3(7 1 + (7 2 ][g(O,s)-0(m c ,S)] , 



where 



c (o) 



3d + C 2 + 3C 3 + C 4 + 3C 5 + C e , 



n 



and the loop function g (m q , s) stands for the loops of quarks with mass m q at the dilepton 
invariant mass s. This function develops absorptive parts for dilepton energies s = 4m 2 : 



g lm«, s 



8, . 8 4 2 . 

— mm a H 1 — v a (2 + y a 

9 ' 27 T 9 V y? 



1 + 



Vq 



1 - 2/ff| 

Z7T ] + ®(y q — 1) 2 arctan 



where m g = m q /mb and y g = 4m 2 /s. In addition to these perturbative contributions, the 



cc loops can excite low-lying charmonium states ip(ls) 
represented by Y LD [p5| : 



t/)(6s) whose contributions are 



Vls s) = — 



x 



a" 



K/Kb £,(0) 



E 



(3C 3 + C A + 3C 5 

TT^r (Yi — ► £ + £~) My t 



Vi=^(u),-,^(fis) ( M Vi - §m b - iM vJ^v t 



where the Fudge factors (see for example 0). 

Let us first study the dependence of the asymmetry parameter olq p on the new Wilson 
coefficients. Note that in further analysis, only short distance contributions are taken 
into account and integration over q 2 is performed in the full physical region 4m 2 < q 2 < 
Omg — my) 2 . We assumed that all new Wilson coefficients Cx are real, i.e., we do not 



m v 

introduce any new phase in addition to the one present in the SM. 

In Figs. (1) and (2), we present the dependence of ag p on the new Wilson coefficients, 
for the B — ► 7r + 7r~e + e~ and B — ► 7t + tt~t + t~ decays, respectively. Here and in all of the 
following figures, zero value of new Wilson coefficients Cx correspond to the SM prediction. 
In the case of B — > 7r + 7r~e + e _ decay the asymmetry parameter ag p is more sensitive to 
C L °l and Crl, while for the B — > h + -k~t + t~ decay it depends strongly on Crl- These 
dependencies can be explained as follows. For the B — > 7r + 7r~e + e _ decay, if the terms 
proportional to electron mass are neglected, it easily be seen from Eq. that 



a 6l 



N2 



- 1 . 



In the limit v — > 1 we get 



No 



+ 



2 (C e 9 ff + C,o + Ct l R ) H - 4C e 7 ff ^H + C RR h ( 
2 (c 9 e// - Q + C L ot L ) H - 4C e 7 ff ^H + C RL h 



(33) 



In the SM in the large dilepton mass region, say about q 2 ~ 5 GeV 2 , Cg + C\o 



Re 



U 9 



c 



10 



9.5. It follows then that the interference terms between Ct/* 



0.4 and 
C10 and 

Cll (Crl) are dominant and hence contributions coming from Cll (Crl) are large. These 
figures illustrates that the contributions of Cll and Crl to ag p is positive for Cll > and 
Crl < 0, and negative for Cll < and Crl > 0. It should be noted that the asymmetry 
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parameter ag p can get only positive or negative values for the case Cll 7^ and Crl 7^ 0, 
while it is always positive for all other choices of the Wilson coefficients, as is the case in 
the SM. For this reason determination of the sign of a$ p can serve as an efficient tool for 
establishing new physics. 

In the B — > 7t + tt~t + t~ process however, the situation is slightly different compared to 
that of the B — > 7r + 7r~e + e~ transition. Largest contribution in this case comes from Crl 
and contributions from all other Wilson coefficients are comparable to one another. This 
observation can be attributed to mass of the r lepton, for which ag p is positive for the 
choice of each individual new Wilson coefficients. 

Depicted in Figs. (3), (4) and Figs. (5), (6) are the dependencies of the asymmetry 
parameters ag and f3$ on the new Wilson coefficients Cx, for the B — > 7r + 7r~e + e~ and B — > 
7r + 7r~T + T~ decays, respectively. Figs. (3) and (5) depict that the asymmetry parameter ag 
for the e + e~ channel depends strongly on the new Wilson coefficient Crl, while it displays 
similar behavior for all new Wilson coefficients for the t + t~ case. We observe that the 
asymmetry parameter [3g depends strongly on Crl in the e + e~ channel and on Crl and 
Crr in the t + t~ channel. It is interesting that fig changes its sign when Crl > 2 in 
the e + e~ channel, while it is negative for all values of C RL or C RR in the t + t~ channel. 
Therefore determination of the sign of f3g is useful in looking for new physics. 

In Figs. (7) and (8) we present the dependence of the asymmetry parameter f3 9 for the 
B — > 7r + 7r~e + e~ and B — > 7t + tt~t + t~ decays, respectively. From these figures one notices 
that this asymmetry parameter is quite sensitive to the variation in Crl and changes its 
sign at Crr > 2 for both channels. Our investigation of the dependence of the asymmetry 
parameter a^ on the new Wilson coefficients shows that for the range —4 < Cx < 4, a v 
varies between —6 x 10~ 3 to 6 x 10~ 3 for the B — > n + 7Y~e + e~ and —5.0 x 10~ 3 to 2.5 x 10~ 3 
for the B — > tt + tt~t + t~ decays, repectively. Therefore detection of the dependence of the 
asymmetry parameter a v on the new Wilsom coefficients is quite hard from experimental 
point of view. 

The asymmetry parameter A\ shows strong dependence on Crl and Cll for the e + e~ 
channel, as depicted in Fig. (9), whose contributions are dominant compared to the other 
Wilson coefficients. For the r + r~ channel contributions of Crl and Crr become dominant, 
as can be seen in Fig. (10). The asymmetry parameter A 2 varies considerably for the e + e~ 
channel, in relation to the variations occurring in Crl, while this behavior is switched to 
the Wilson coefficient Clrrl for the r + r~ case which are presented in Figs. (11) and (12). 

Presented in Figs. (13) and (14) are the dependence of the asymmetry parameter 
on new Wilson coefficients for the e + e~7r + 7r~ and t + t~7t + tt~ decays, respectively. In both 
cases the asymmetry parameter A v shows strong dependence on Crl- 

Finally, in Figs. (15) and (16) we present the dependence of the averaged (i.e., integrated 
over q 2 in the full physical region) CP asymmetry on the Wilson coefficients for the B — > 
7r + 7r~e + e~ and B — > tt + tt~t + t~ decays, respectively. We observe that (Acp) is strongly 
dependent on Crl and Crr for the e + e~ channel, while this dependence is switched to Cll 
for the r + r~ channel. One can easily read from this figures that (Aqp) > —2.0 x 10~ 2 in 
the e + e~ channel, for the values Crr < —2 and |(Acp)| > 1.7 x 10 -2 when Cll < —2. 

As the final concluding remark, we presented in this work the model independent anal- 
ysis of the exclusive B — > (£ = e, r) decay is presented. In particular, the 
sensitivity to the new Wilson coefficients of the experimentally measurable asymmetries 
and CP violating asymmetry are systematically analyzed. The main result of the present 
study is that different asymmetry parameters show strong dependence on different new 



13 



Wilson coefficients. Therefore a combined analysis of the different asymmetries and CP 
violating asymmetry can give unambiguous information about the existence of new physics 
beyond the SM and especially about various new Wilson coefficients. 
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Figure captions 

Fig. 1 The dependence of the asymmetry parameter a$ p on the new Wilson coefficients 
for the B — > 7r + 7r~e + e~ decay. 

Fig. 2 The same as in Fig. (1), but for the B — > it + it~t + t~ decay. 

Fig. 3 The same as in Fig. (1), but for the asymmetry parameter a$. 

Fig. 4 The same as in Fig. (1), but for the asymmetry parameter (3 g . 

Fig. 5 The same as in Fig. (3), but for the B — > it + it~t + t~ decay. 

Fig. 6 The same as in Fig. (4), but for the B — > it + it~t + t~ decay. 

Fig. 7 The same as in Fig. (1), but for the asymmetry parameter (3 V . 

Fig. 8 The same as in Fig. (7), but for the B — > it + it~t + t~ decay. 

Fig. 9 The same as in Fig. (1), but for the asymmetry parameter A\. 

Fig. 10 The same as in Fig. (9), but for the B — > it + it~t + t~ decay. 

Fig. 11 The same as in Fig. (9), but for the asymmetry parameter A 2 . 

Fig. 12 The same as in Fig. (11), but for the B — > 7r + 7r~T + T~ decay. 

Fig. 13 The same as in Fig. (9), but for the asymmetry parameter A v . 

Fig. 14 The same as in Fig. (13), but for the B — > 7r + 7r~T + T~ decay. 

Fig. 15 The dependence of the averaged CP asymmetry on the new Wilson coefficients 
for the B — > 7r + 7r~e + e~ decay. 

Fig. 16 The same as in Fig. (15), but for the B — > 7r + 7r~T + T~ decay. 
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